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ABSTRACT 


Although simple radial equilibrium turbine design 
has long been a useful design method for axial flow gas 
turbine blading, more demanding requirements for high 
performance machinery make such an approximation less 
justified. More exact analytical methods do exist, but 
require many times more effort, even if a computer is 
utilized. 


A relatively simple model for improving a simple 
radial equilibrium preliminary design is proposed in 
this thesis. The axial velocity ‘distribution, determined 
from the input SRE design, is modified to take into 
account the results of flow studies in flared annuli 
and of actuator disk theory. Using the corrected axial 
velocity distribution, new velocity triangles can be 
calculated. 


A FORTRAN IV computer program to implement the 
proposed model was written and is presented. Its 
potential is demonstrated by application to. three 
example simple radial equilibrium designs.:. The results 
Show practical differences exist between a Simple radial 
equilibrium design and the deSign corrected by the 
proposed method. 
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NOTATION 


filenareayeeet— 


SCOcrmpelenG In Couations 2.24, 2.25, 2.27, 
and 2.28 


annulus area, ft? 


SOerrrertenemim Equatzons 2:24, 2.25, 2.27, 
and 2.28 


specific heat at constant pressure, BTU/1bm/°R 
specific heat at constant volume, BTU/1lbm/°R 
Ben? f eelomybt/sec- 


specific enthalpy, BUT/lbm; actuator disk 
height, ft 


Woe it 1 bf 7B0U 

Mach number 

mass flow rate, lbm/sec 

turbine rpm, revolutions per minute 
pressure, pt /it- 
GetewOmmiieatetramsi:er to fluid, BTU/sec 
radial distance from machine axis, ft 
Srpecitgesoas Constant, ft ibf/lbm/°R 
specific entropy, BTU/1bm/°R 
temperature, °R 

blade velocity, ft/sec 
Pluscevelocity, ft/sec 

rate of output shaft work, BTU/sec 
relative velocity, ft/sec 


axial distance from an actuator disk, ft 





fo absolute angle measured from axis, degrees 


Aa Stator turning angle, degrees 

B relative angle measured from axis, degrees 
AB rotor turning angle, degrees 

yi specific heat ratio 

0 Consity welbm/ ft - 

d flare angle, degrees 

w) angular velocity, radians/sec 
Subscripts 

h hub, Comal 10n 

i reference condition 

O EOUQmmCSweagmMattom) COnd1 tion 

ag radial component 

rel relative 

= ti pecenali lon 

x axlal component 

0 tangential component 

1 Geueteiinlec, Plane one condition 
2 BOuOnpweouEele:; plane two condition 
3 plane three condition 

+00 condition far downstream 

—o0 condition far upstream 


NOTE: Notation not appearing here may be computer 


program variables listed in Appendix A. 


10 





eet Nee ODUCTION 


1.1 General 

Energy is a topic of increased interest in current 
times. Searches for new sources of energy and new means 
of converting energy into useful Rogen have been under- 
taken, along with efforts to conserve the resources 
available. While exploiting solar energy, tapping 
geothermal power, and drilling for oil off-shore receive 
mucma public interest, progress in improving existing 
conventional forms of energy transfer and conversion must 
also be made. The great importance of energy utilization 
has made it worthwhile for designers of energy related 
systems to pay even more attention to achieving the best 
possible performance. Improving the efficiencies of the 
various hardware components which carry out the heat 
transfers and energy conversion functions of. the thermal 
cycle is one desirable goal. To econo sneen & eS. 
useful to make design methods more accurate and more 


easily applied to real problems. 


1.2 Axial Flow Turbines 

This thesis deals with one commonly used component 
for converting energy into useful work--the axial flow 
turbine. Turbines convert the kinetic energy which has 


been provided the working fluid into mechanical work 


Ju 





Perrermed upon a rotating shaft. Axial flow turbines are 
so named because the flow of the working fluid through 
the machine is predominantly in a direction parallel to 
the machine's axis of rotation. Such turbines can 
Operate with a variety of working fluids, and in either 
open or closed thermal cycles. The most common con- 
figurations are aS a steam turbine in the Rankine Cycle 
(a closed cycle) or aS a gas turbine in the Brayton 
evelert{an Open Om closed Gycle). Historically, the 

steam turbine found earlier application than the gas 
turbine, which recently has gained in importance. The 
development of the gas turbine depended upon advances in 
high temperature material properties and was also closely 


tied to the successful operation of axial compressors. 


1.3 Aerodynamic Design of Axial Flow Turbines 
1.3.1 General 

One aspect of the design, which procedes similarly 
for axial turbines of all kinds, is the aerodynamic 
design of the turbine blading. The efficiency of the 
turbine depends upon the smooth transfer of kinetic 
energy from the fluid flow, with as little energy wasted 
aS possible. The shapes and angles of the turbine 
blading must be appropriate to the pattern of the flow. 
The aerodynamic design of the turbine determines what 


these flows and shapes should be. 
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In general, the aerodynamic deSigner of an axial 
turbine follows the steps listed below: 

(1) for an assumed operating point and 
efficiency, carry out the thermodynamic analysis of the 
thermal cycle to determine the mass flow, and inlet and 
ouimmet States for the turbine: 

(2) postulate a working fluid flow which 


satisfies the laws of motion and produces the required 


work; 

(3) design the shape of the turbine blading 
appropriately; 

(4) calculate the efficiency and evaluate the 
design; 


(5) repeat all of the above steps as necessary. 

ieee ee loweso lutions 

As is true of most design efforts, there are many 
levels of sophistication and detail which the designer 
may seek, depending upon the stage of the design ang / OG 
the cost and criticality Of ee Product. In the aero- 
dynamic design of a turbine, the designer must decide 
which level is best in each case. For, although the 
equations of motion and continuity for a fluid are 
readily written, they cannot be solved in practice unless 
simplifying assumptions are made. Fortunately, such 


Simplified forms have proved adequate in many applications. 
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1.3.3 Two-Dimensional Flow 

Pemenc no -Etomratlo Of a turbine 1s nearly 1.0, 
no large radial component of fluid velocity is expected, 
and flow conditions will be Meeronmnate ty the same from 
root to tip. The flow can then be analyzed in two 
dimensions. This results in blades of constant cross- 
section, which have been extensively studied experimentally 
in the form of cascades. 
1.3.4 Three-Dimensional Flow 

In some cases the designer faces requirements on 
blade tip speed, pressure ratio, and mass flow which 
could result in low hub-tip ratios (i.e., relatively 
long blades) and considerable turbine flare. These 
conditions reduce the accuracy of simple two-dimensional 
design methods and force the consideration of three- 
dimensional flow. The general solution is impractical, 
but simplified methods have been developed and used with 


Success. Several of these approaches are mentioned below. 


1.4 Approaches to Three-Dimensional Design 

The Bae assumption made in virtually all cases is 
to assume axial symmetry. The number of blades is 
(Var tolenelywlarce that Circumferential conditions at a 
station can be represented by average vatues.— 
ede weeaaeed Lb hou do bra um 

Assuming axisymmetric flow, the differential 


equation of motion in the radial direction for the fluid 
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can be put in a form (see paragraph 2.4.1) which may be 
numerically integrated by a computer. The procedure in- 
volves complicated iterations, however, and elaborate 
measures must be taken to ensure the proper convergence 
of the streamline locations. Carmichael has discussed 
the procedure and presented block diagrams of possible 
computer programs. ~ This method allows for fluid 
velocities in the radial Bli@ection whiten cresudser from 
imbalances in the static pressure and centrifugal 
acceleration forces on the fluid, but its complexity 
1s formidable. 
im 2 simple Radial Equilibrium 

If the radial components of the fluid velocity are 
neglected, the so-called streamline curvature term in the 
differential equation drops out. The differential 
equation of motion in the radial direction can then be 
written involving only the axial and tangential velocities. 
This can be easily integrated for axial velocity when 
given a reasonable distribution for the tangential 
velocity (see.paragraph 2.4.3). Simple pairs of functions 
for axial and tangential velocity, so determined, have 
been widely used in three-dimensional design. 
1.4.3 Actuator Disk Theory 

Simple radial equilibrium analysis depended upon the 
radial fluid velocity being zero before and after each 
blade row. In most cases, changes in the radial distri- 


Purton or Genstey result am radial shifts in the 


iS 





“wednlenes,  Thererore, a non-zero radial velocity 
component must exist. 

An alternative to the simple radial equilibrium 
analysis is the actuator disk formulation. Each blade 
row 1s replaced by a narrow “actuator disk" which provides 
the appropriate anatan change in fluid tangential velocity. 
This device permits an approximate solution to the 
equations of motion, allowing finite radial velocities, 
but still is restricted to annular flow between walls 


parallel to the machine axis. (See paragraph 2.5) 


1.5 Fluid Flow in Flared Annuli 

A frequent deSign practice is to keep the axial 
velocity approximately constant throughout the machine, 
compensating for changes in fluid density by increasing 
the annulus area. This area change results in "flare"; 
Bidets, ene wnubeanad shroud Surfaces are not parallel to 
the machine axis, but are sloped. 

This flare could be expected to cause a general 
change in the axial velocity distribution within the 
annulus, with significant radial velocity components 
appearing near the walls. Carmichael and Pai“ have 
Suggested a relation which can be used to calculate the 


change in axial velocity in flared annuli. 


1.6 Summary 
The preceding discussion was intended to demonstrate 


a need for a relatively simple method to predict the 
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fluid flow distribution in the presence of flare, improving 
upon the accuracy of the simple radial equilibrium model. 
Such a method would hopefully be useable in the preliminary 
design phase with little increased effort, as compared to 
the more complicated streamline curvature approaches. 

A method is proposed in this thesis for "correcting" 

a preliminary simple radial equilibrium turbine design 

for flare and streamline curvature effects. This is 
accomplished by modifying the axial velocity radial 
variation, as determined from the simple radial equilibrium 
equations. A correction for turbine flare is applied 

using the findings of Pai; next, the axial velocity is 
further modified at the blade leading and trailing edges 

by the approximate solutions to the actuator disk analysis. 
As the computations are Still iterative in nature, they 

are very suitable for computerization. 

A computer program to modify a simple radial equili- 
brium turbine design, according to the model described 
above, is presented in this thesis. While not theoreti- 
cally deep-seated, it is hoped the methodology may be 


helpful in some turbine design applications. 


uae 





2. BACKGROUND 


The purpose of this chapter is to provide the back- 
ground for the methodology employed in the computer model 
to be discussed in Chapter 3. The emphasis is placed on 
those assumptions and relations which directly apply to 
the program model, although some additional development 
1s included for perspective and to point to possible 


areas of application and further work. 


2.1 Fluid Mechanics 

Axial turbine flow is frequently modeled as incom- 
pressible, at least in the radial direction, so that the 
equations of motion for an incompressible fluid apply. 
It is also known that viscous effects can usually be 
neglected outside the boundary layer.? Pp cy Pica eeark 
coordinates, the equations of motion for an ideal, incom- 


pressible fluid can be written 2ao 





G 
ig 0° xr 1G GC @ eps 
mr arco x Ox eal) ee eet) (2o8, 
ooo oe, fe feta Se oP (232) 
ie Be ec 0 X aX i or 08 : 
Wee eo Vx _ _ Jo ap ec) 
tor r o6 X OX Of a es 
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where Ve Vos and Ms refer to components of fluid velocity 
in the radial, tangential, and axial directions, respect- 
ivelywest lS the Static pressure. The continuity 


equation is 


oV. Bs OV avi 
ae Go ay (Boe) 


a Le 


2.2 Thermodynamics 
Zo ee NOLK 

From the first law of thermodynamics, the steady 
flow energy equation, describing the change in state of 
a fluid system in flowing through a control volume, can 


Bee written ae 


Me 7 v* 


Q- We=m{[(h + 7907 out (h + 2g9 


). J], (225) 


1g 
where potential terms have been ignored. QO is the rate 
of heat transfer to the fluid, W is the rate of output 
shaft work, m is the mass flow, V is the fluid velocity, 
and h is the enthalpy. A useful definition is to make 


Stagnation enthalpy, Roe 


2 
Wo eee 2 


. 253 (2.6) 


Using equation (2.6), and assuming that the heat transfer, 
Oo, from an axial turbine is small (compared tO W), 


eguation (2.5) becomes 


le, 





(25 oe) 


where the 2 


ke 


ana I refer to the outlet and inlet states 
respectively. 


Following reference 1, the rate of output shaft work, 


W, is also equal to the product of the angular velocity, 


Ww, Of the blades and the torque caused by the moments of 


the external reaction forces. The torque is given by the 


change in angular momentum of the fluid, or 
Torque = m(rjVo5 - 


r,Vg))- (2 ey) 


Since 


ste 
l 


LOGeguc = W& 


and the blade velocity, u, is just wr, then 


Wee eo. 1 Ol 
Ee aw —- ee a (2.9) 
ry ze; a | 

Mm O 


In axial machines. it is a good assumption to make U5=U, 


i@~e preliminary design. 


2.2.2 Gas Laws 


The equation of state for a perfect gas can be 


Written 
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P= ORI, (27r OD) 


where R is the specific gas constant. The following 
relations between the specific heats are the standard 


ones. 


a 
oe a 
V 
C., GS Ss 
Oe as 
C,, (<a) 7 (2.1) 


f.2.o) The Second Law 
From the second law of thermodynamics, uSing Gibbs 


Equation and the perfect gas law, the specific entropy, 


Ss, can be written agtye 
fda & dh 2 =. aie (2.12) 
OT : 


2.3 Compressible Flow Relations 

In all but very low speed turbomachines, in which 
velocity, pressure, and temperature changes are very 
small, the effects of compressibility of the working 
fluid become important. Consideration of three-dimen- 
Sional compressible flow is extremely complex. Fortunately, 


the changes in fluid density are often not great in the 


au 








radial direction at a station, and need to be considered 
only in the axial flow direction. 

The one-dimensional compressible flow relations for 
an ideal gas appear in many representations and are 
conveniently used. Horlock> has a discussion of their 
derivation and summarizes the results. The following 
forms will be useful in discussion later. The meanings 
of the symbols can be found in the notation list. From 
equation (2.6), since, for a perfect gas, h = Sa 

T, = P+ oe (Zale) 

O Pp 
Homngq (2.11) and dividing by Tor equation (2.13) can be 


GewVerttcen as 


= 2 
z=1- (+) (2.14) 
O 2gR To 
The Mach number, M, 1s defined by 
Mies a (2.15) 
Se 
g.YRT 
O 
PaGmmeouatironms (2,14) and (2.15), 
ly cpeey 7 
eee eee et (2.16) 
M v2 Z 
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Srewpressure-temperature relation for an isentropic 
process can be shown to apply in the compressible case, 
and gives 


ae 


p = (=) Y*, opale)) 
O O 


From continuity, the flow area, A, is 


ni 
Pi eomeea ie Vewdetemuincad trom equations (2.14) and (2.17) 


above. 


2.4 Simple Radial Equilibrium 
2.4.1 General 

As indicated in the introduction, when the hub-tip 
ratio of a turbine stage is further and further from the 
Volcom lam loc mmmmelatavely Jonger and longer blades), 
the two-dimensional flow design methods become less and 
less adequate. In some marginal cases, three-dimensional 
details of the design are ignored in favor of other 
considerations, for example, maintaining constant blade 
sections. But modern performance requirements make such 
practices less attractive. The general three-dimensional 
flow equations are difficult to solve, however, and 
Simplifying assumptions are frequently made. The most 


general such assumption is to neglect any variation of 


Zs 





Ce eee 





properties in the circumferential direction. In the 
feral Crecerion, where tluid property changes may be 
large, the one-dimensSional compressible flow relations 
are assumed to describe the flow. Finally, in the radial 
direction the equation of motion in a simplified form-- 
termed the simple radial equilibrium (SRE) equation--is 
often applied. A discussion of the development and use 
of the simple radial equilibrium equation follows. 
2.4.2 The Simple Radial Equilibrium Equation 

Assuming axial symmetry, the equation of motion in 


the radial direction, equation (2.1), can be written as 
Vom t+ VO Oe Or EE rl (27.09) 
De 


a + Ve and by differen- 


tiating with respect to r, equation (2.6) becomes 


From the identity v2 = Vee + V 


oV oV oV 


a eee, ae _ = ae __* 
ox tC« el eee Go sar | x or! Ga 


and equation (2.13) can be rewritten in the form 


Q? 
1e | 
Q? 
” 


i oh 
—— eee 2 ee -_ ——, 2 e Z 
oJ or z oc ox ( a 


Substitution of equations (2.20) and (2.21) into equation 


(2.19) gives 


as O ig x 8 8 (2.22) 


aa or x OX x OL 6 or r 
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MeehtS point dS/dr and dh fox are often assumed to be 
Zero. ACcCOrding to Horlock® these two assumptions, 
which prove to be nearly true in practice, hold for the 
flow between blade rows of a reversible turbine in which 
delivered work and total pressure drop across the blade 
are the same at all radii. The tern, Vi. OV 79x, 1S 
called the streamline curvature term and is the cause of 
Piemolulcwety,—temaLoned in the introduction, in 
obtaining a solution to equation (2.22).> In addr ti1e68n 
to assuming ds/odr and dh fox are equal to zero, if it is 
assumed that ve equals zero at the leading and trailing 
edges of a blade row, the streamline curvature term drops 


out and the partial derivatives can be replaced by 


ordinary derivatives, giving 


dv dv ae 
ee 8 


ear e) dr r ae 


Equation (2.23) is called the simple radial equilibrium 
equation. 

A more direct derivation of the simple radial 
equilibrium equation, which may be more Intute ve,- Canmee 
developed from consideration of an infinitesimal volume 
of fluid in equilibrium with the external forces acting 
on Ae The details are shown in Figure 2-1, on the 
next page. The substitution for the g/p (dP/dr ) term 


procedes as above. 
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Figure 2-1 
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In summary, the simple radial equilibrium equation 
Toman for che following conditions: 
(1) ideal, incompressible fluid flow 
(2) axial symmetry 
(3) perfect gas law holds 
(4) Ve = 0 before and after blade rows 
(5) dh /dr = 0 
(6) ads/dr = 0 
Peee oe cOlutions to the Simple Radial Equilibrium Equation 
From equation (2.9) it can be seen that the work 
performed by an axial turbine is related to the tangen- 


tial velocities before and after the rotor, and V 


Yoq 62° 
respectively. A design problem for an axial turbine is 
to propose a reasonable tangential velocity distribution 
which will give the desired work, and then by means of 
the simple radial equilibrium equation find the axial 
velocity, Vier so that the blade shapes can be determined. 


For simple functions of V, = Vo (x), equation (2.23) is 


8 
eaSily integrated directly for cee = V(r). 

It would appear that independent functions for Vnq 
and Vo5 could be chosen and corresponding functions 
Vi7 (Y) and Vig (4) calculated from the simple radial 
equilibrium equation. However, the Vy's are related by 
the work relation, equation (2.9), and by the assumption, 
used in deriving the simple radial equilibrium equation, 


that the stagnation enthalpy, Aas after each blade row, 


does not vary with radius. 
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A familiar, useful set of tangential velocities 


given by Carmichael and Lewis, and described by Horlock®, 
are the following: 
Ss n b 
oles oe” (2.24) 
Le ia ys ; 
Vo2 = ar as (2.259) 


where a, n, and b (b > Q) may be chosen by the designers. 


Maen, fLrom equation (2.9) 


ve Vaiba\\| 
02 60g J 





(V cai (V 


01 ~ Yeo): 


and using equations (2.24) and (2.25) 


W _ 2mrN ,2b, _ 4mNb 
= 0G a ens, ae oe 
m O O 


Equation (2.26) shows the work done (hence hy) 1s 
independent of radius and has the proper sign for a 
Bur bane . 

A modification of the set of tangential velocities 
(2.24) and (2.25) is considered in this thesis. Suppose 
a set of tangential velocities of the following form is 


given: 


28 
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_ 2 a 

Vor = aor + ayr + ay + =a (Zoe) ) 
2 Do 
Vg = a,x a a,x ta, + aor (2c ee) 
Here also the work is independent of radius, for 

i 2a ee (2229) 
° 60g_J 1 2° ' 
m Oo 


Here by must be greater than b. if work is to be done by 


2 


the turbine. The constant by is virtually always positive, 


and b. 1s usually negative or zero. The advantage of 


allowing by and bo to be different is to be able to 


provide for special cases such as ‘free vortex, axial 
outflow (25, Ayr ays ba = 0: by 0). 

Substituting an expression of the form of (2.27) 
into the simple radial equilibrium equation and integrating 
from reference point eer = V(r; ) EOr pei nt ae = Viale the 


following expression for axial velocity is obtained (the 


details are in Appendix D): 


10 as 8 
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Equations (2.27) and (2.30) for tangential velocity and 
axial velocity form a set which satisfy the simple radial 
Poueelomilun equattommi2.23). It is important to note 
that Many factors and various criteria for good practice 
influence the detailed choice of the constants in the 
expressions for Vo- 
Two well known special cases of equations (2.27) 
and (2.28) are the constant reaction Vo qistribueicon and 
the exponential Vo distribution. These tangential 
velocity distributions are presented in Avpendix D 


together with the corresponding function for axial 


velocity. 


2.5 Actuator Disk Theory 
2.5.1 General 

It is recalled that the basis for the simple radial 
equilibrium design method is that the radial component of 
the fluid velocity be zero at the leading and trailing 
edges of a blade row. This assumption results in stream- 
mines sok the form 2ailustrated in Figure 2-2(a). A more 
general method of analysis, which takes some account of 
radial shifts in the streamlines between blades, is the 
actuator disk model. As described by Hover a rotor 
or stator row is replaced, conceptually, at approximately 
midchord, by a thin disk across which the tangential 
velocity change takes place. Such a device allows the 
equations of motion, including radial velocity terms, to 


be solved analytically, although the formulation is still 


complex. 
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moe ADDEOximate Solution 

With some simplifications, an approximation to the 
complete actuator disk results can be obtained. Radial 
WelOecities are assumed to be zero at the walls and small 
elsewhere, and Simple radial equilibrium is assumed to 
Pemvalid a surficient distance from the disk. According 
to this model, the streamlines will be similar to those 
shown in Figure 2-2(b). The approximate solutions for 
the axial velocity near a Single disk, in terms of the 
axial velocity on the same streamline far up and down- 


stream, are: 


a (val. 
v(x) =(V_) _ + a3 = exp (—*), (223%) 
x x“ 00 D h 
Hestieamn OL stne arsk (x < 0), and 
(V_) = (ae) 
= S22 eee See ales 
V(x) = Ore ae : exp ( ell C2432 } 
Ao umemaean Ormenewamsk (x > 0). The notation Ee ae 


and (i lee refer to the axial velocities given by simple 
radial equilibrium a long way downstream and upstream, 
respectively; h is the height of the actuator disk. 

This analysis, like the simple radial equilibrium 
analysis described in the previous section, treats ideal, 


incOnemesstole ntlow in cylindrical annuli. 
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For real turbines, where many rows of blades (disks) 
may influence the axial velocity distribution, the effects 
of each disk may be considered separately, then super- 
imposed to find the real axial velocity.* This procedure 
is used in Appendix E to find the axial velocity distri- 
bution at a point influenced by two nearest adjacent 


disks. 


Zao Flow an Flared Annuli 
2.6.1 General 

Mmad toto LOenproving the methods for cylindrical 
flow, efforts to understand the fluid flow in ducts of 
changing annular cross-section have been made. Such 
sections are commonly used in practice. 

Horlock, in reference 6, describes some theoretical 
models, proposed by various workers, for the determination 
of axial velocity. Analyses by Lewis and Horlock for 
incompressible flow in conical ducts, by Wu for compressible 
flow in flared annuli, and a simplification of Wu's analysis 
by Walker are discussed. The drawback of these analytical 
methods is their complexity. 

Experimental work by Carmichael and ea has 
suggested a simple expression relating the axial velocity 
of flow in a flared duct to the angle the streamline 
makes with the axis. Their work, Gescribed below, forms 
the basis of the flare correction to axial velocity 


given in this thesis. 


33 





2.6.2 Experiments Representing Fluid Flow in Flared 
Turbines* 

Experiments were conducted using a tilted-floor 
pilecmerolvyeEic tank to measure flows within a flared 
annulus. An electrolytic tank provides an analogue to 
ideal incompressible fluid flow. Electrical potential 
measured in the tank 1S analogous to fluid potential and 
the electrical potential gradient is analogous to the 
fluid velocity. A sector of an axisymmetric model can 
be represented by tilting the floor of the tank. 

Although flow in a real turbine is compressible, 
the radial variation in density at any axial station is 
relatively small. Consequently, it can be assumed that 
in usual geometries the seogenllane positions for the 
real, compressible fluid are very nearly the same as for 
the ideal, incompressible case. 

One aspect of the fluid flow measured was the 
general change in axial velocity distributions in the 
flared annulus. Three models were studied, with 
different distributions of flare between the hub and 
shroud walls; each model had a hub-tip ratio of 0.9. 
fjemolctributicon was expected to be a function of location 
within the annulus, described by an angle, ¢, defined as 
shown in Figure 2-3. The flow conditions were measured 
at an axial position, free of discontinuities, such as 
section A-A' in Figure 2-4. The axial velocities were 


found to have the distribution shown in Figure 2-5, 
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approaching the flow from a source. The equation, 
Peeoring tO pPetential flow theory, for the axial 


Melocity drstribution £or a two-dimensional source is 


V 


x ae 2 
Vi. ($=0) me COLE E (255 8) 81) 


Pai suggests that comparison of this equation, also 
plotted in Figure 2-5, with the experimental results, 
Hicdchertes that for angles of ¢ < + 15°, equation (2.33) 
1S a good representation of the axial flow. Even at 
larger angles, the agreement may be good enough to use 


equation (2.33) for some design applications. 


2./ Axial Turbine Flow Geometry 
Zeya, Velocity Tans 

It should be restated that the purpose of deter- 
mining the details of the fluid flow is to be able to 
design the blades appropriately. The relevant fluid 
velocities and flow angles are conveniently displayed 
in the form of vector diagrams, often called velocity 
triangles. Velocity triangles for a single stage axial 
turbine with pure axial inflow are shown in Figure 2-6, 
where any radial velocity components have been neglected. 
iiomealeulation Of Such velocity triangles for each 
leading and trailing edge and at several radii must be 
accomplished before the detailed design of the turbine 


blading can continue. 


a) 








Stator anlet Vv 






ROLOr anleL 


Rotor outlet 


Remarks: 


1. the signs of velocity corponents are positive in the 
direction of U or V, 


2. angles are positive in a counterclockwise direction 
from Ve 


3. as drawm above-- 
Ve1s X1, A, are positive 


Veo, Xo, Po are negative 


Figure 2-6 
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2.7.2 Annulus Area 

Suppose a fluid is flowing uniformly through a small 
area at a given rate. If the area encompassing the flow 
is perpendicular to the streamlines, it is termed the 
Muorvmatea, A, and 1S given by equation (2.18). I£ the 
area is in a surface whose perpéndicular makes an angle 
@ to the streamlines, as in Figure 2-7, then the area 


required to encompass the same flow is 


oem = A/cos a. (2.23549) 


In the case of an axial flow turbine, see represents 
the annulus flow area perpendicular to the machine axis, 
where @ is the same a as in Figure 2-6. From Figure 2-7 


it can be seen that 


= Vv: 
cos a aA 


whence equation (2.18) becomes 


A == ean (2.2354 
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Fisure 2-7 
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3. MODEL DESCRIPTION AND IMPLEMENTATION 


3.1 General 

In the introductory chapter the need for a simple 
methodology to generate a preliminary axial turbine design, 
taking account of turbine flare and radial streamline 
shifts, was discussed. The simple radial equilibrium 
design method was presented in Chapter Two, together 
with the results of actuator disk theory and axial flow 
distribution studies in flared annuli. In this thesis 
a procedure is proposed which accepts a Simple radial 
equilibrium preliminary design as a basis, and then 
modifies the axial velocity distribution in accordance 
with the actuator disk and flare corrections mentioned 
above. These modifications to the simple radial equili- 
brium flow solution result in slightly different blade 
shapes, which should be more commensurate with the real 
flow. 

This chapter will present the model for modifying a 
Simple radial equilibrium preliminary design, and describe 


a computer program which implements the model. 


bezertiode! Description 
To begin, it is useful to consider a sketch of the 


geometry and then a bare outline of the problem. 
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Figure 3-1 represents a Single stage axial flow turbine. 
fe-meeacor inlet, rotor inlet, and rotor outlet planes 
are labeled 1, 2, and 3 respectively. Other features 
will be discussed as they become important. 
The following are steps in the procedure which may 

be identified in order to guide the discussion: 

CE Input mem scimple radial equili bream 
preliminary design 

(2) at each of the three principal planes of the 
turbine stage, calculate a set of streamlines which is 
consistent with certain given conditions and the 
modified axial velocity 

(3) at every streamline in each principal 
plane, calculate velocity triangle data. 
The rest of this eee will discuss the above steps in 
general terms. The details of the procedure are best 
wum@erssceood in the context of the computer program des- 
Satecien, beqinning in paragraph 3.3. 
3.2.1 The Input Design 

The purpose of the procedure discussed in this 

tmesis 1S to make slight “corrections” in an existing 
Sipe radial equilibrium design in order to refine the 
detailed blade shapes. As a result of carrying out a 
preliminary design for a single stage, many parameters 
will have been established. The working fluid, total 
temperatures and pressures, hub and tip radii, machine 


rpm, mass flow, and blade widths are known; also, in the 
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@2-c of Simple radial equilibrium, an appropriate pair of 
functions for axial and tangential velocity at each aiawe 
will have been picked. These quantities are indicated in 
mere 3-1. 

For present purposes, it is assumed that it is 
desired to retain the input annular geometry (hub and 
mueesadii, and blade widths) in the modified case, 
although other choices are possible. For example, the 
average axial velocity could be fixed, and the tip radius 
allowed to change to accommodate the flow. 

Horn simplicity, the input preliminary design is 
restricted to having hub and shroud surfaces of uniform 
Slope over the whole stage, aS in Figure 3-l. Also, the 
rotor tip clearance and stator-rotor gap are assumed to 
be zero. The model could be applied to more general 
cases with relatively minor modifications to the computer 
program. The turbine flare may be distributed between 
the hub surface and shroud surface in any proportion. 

It should be noted that the procedure described 
here will, in no sense, improve a bad preliminary deSign. 
3.2.2 Calculating the Streamlines 

To fully determine the blade shape, fluid flow 
Moreormation, including axial velocity, must be known at 
several radii. If the design were merely Simple radial 
equilibrium, the velocity triangles could be independently 
found at any convenient radius, Since ie and Vo (and U) 


are known functions of r alone. However, in the present 


44 








@---,, the actuator disk correction to the axial velocity 
at a location on one plane is dependent upon the axial 
velocity along the same streamline at stations up and 
downstream. For this reason a set of streamlines passing 
ali three principal planes must be calculated. 

To find the final set of streamline positions, three 
different iterative loops are involved. They are: 

(1) at a given plane and streamline, iteration 
is required to obtain the next (radially outward) 
streamline position; 

(2) at a given plane, iteration is required to 
obtain an average axial velocity which will result ina 
match between the fixed mass flow rate and annulus area; 

(3) since a change in the axial velocity ata 
Streamline in one plane affects the axial velocity in 
another plane, the streamlines are recalculated at each 
plane in succession until the changes from one iteration 
to the next become inSignificant. 

The modified axial velocity, used in the streamline 
eeveulatcions, 1S found utilizing the approximate actuator 
disk results described in paragraph 2.5 and Appendix E, 
and the flared annulus axial flow relation described in 
paragraph 2.6. The only changes to the simple radial 
equilibrium preliminary design are a result of these two 
meomrcctions' to the axial velocity, which impact on the 


velocity triangles. The changes represent a departure 
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fomesimple radial equilibrium. The justification for 
this procedure is that the adjustments are small and are 
improvements in the direction of the real solution. 
3.2.3 Velocity Triangle Data 

After the positions of streamlines have converged, 
all relevant parameters in the velocity triangles can be 
ea Cusated from =thespraimary variables U, corrected Ves 
and Vor eile: shgyenvis eae RE. Other useful design data, 
fon example, stator and rotor fluid turning angles and 
W/W, ratio, are easily derived from the basic velocity 
triangle data. 

For purposes of comparison, sets of velocity tri- 
angles can be computed at the same streamline-plane 
Pecations, but using the unmodified simple radial 


equilibrium axial velocity. 


3.3 Computer Implementation of the Model 
3.3.1 General 

The iterative nature of the detailed calculations 
virtually precludes hand calculations of the streamlines 
as a practical matter. A computer program to accomplish 
the calculations has been written in FORTRAN IV language 
making use of the model easy and inexpensive. 

The goal of the work described in this thesis was 
womcdevelop the model and to program it, so that it 
could be applied in future design studies. Since the 


model itself is experimental, the program is not merely 
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the mechanization of standard, accepted practice. For 
this reason, it was thought advisable to include here a 
sufficiently detailed description of the computer program 
so that the model used is clear and so that changes can 
be made easily in the future if desired. References 7 
and 8 were useful aS programming language guides. 

The programming philosophy was to simplify the 
understanding and use of the program and to provide 
flexibility. The program was seqmented into subprograms, 
Sfemenaving a Specific function. Identical variable 
labels have the same meaning in all segments of the 
program in order to simplify understanding it. The use 
of COMMON reduces the storage demands of the program. 
Extensive execution information, such as the number of 
loops performed and intermediate values can be printed 
out at the user's option. Diverging iterative loops are 
terminated by internal testing. 

The following paragraph lists the significant steps 
in the execution of the program to modify a simple radial 
equilibrium design in accordance with the model. Sub- 
sequent paragraphs will discuss each step in detail. 
Bee-2 Procedural Steps 

The £Ollevwme list contains a list of all significant 
steps in the program. Some details, such as initializa- 
tion of variables and testing for loop convergence, are 
omitted and are self-explanatory or can be conveniently 


discussed separately. Free use will be made of the 
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notation employed in the computer program itself. The 


list 1s intended to serve as a guide for the discussion; 


full appreciation of the interrelationships can be gained 


only by studying the program flow charts. 


(1) Read and process input data 


(2) At the first plane calculate a set of 


(3) 
(4) 


Streamlines satisfying continuity using 

the simple radial equilibrium axial 

velocity corrected for flare 

a) the first streamline follows the hub 

b) guess at the location of the next 
streamline and iterate until continuity 
1s satisfied between two streamlines 

c) find locations of succeeding streamlines 
in turn until the last one has been 
calculated 

d) compare the radial location of the last 
streamline to the desired tip diameter 
at the plane under consideration 

e) if a discrepancy exists, modify the 
taverageuwmaxlal velocity at that plane 
and repeat from a) above 

Repeat (2) for planes 2 and 3 

After a set of streamlines is calculated 

for all three planes, recalculate them all 

as before but using the actuator disk 


correction on the axial velocity as well 
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(5) Test to see if the streamline positions 
for plane three have changed much; if they 
have, go back to step (4) and recalculate 
the streamline positions for all three 
planes 
(6) Calculate and print velocity triangle data 
for each streamline position at each plane, 
uSing the modified axial velocity 
(7) Repeat (6), but using the simple radial 
eGuiliparum axial velocity. 
Sees Overall Program 
The overall program is represented by the flow 
eiemet in Figure 3-2, on the next page. The circled 
numbers provide points of reference for more specific 
meewecharts in later Eigures. The initialization, 
Poememeniting, and printing of £1, L2, and L3 are included 
in the flow charts so that the user can readily interpret 
Certain diagnostic output to be described later. 
3.3.4 Reading and Processing Input Data 
Memmi ce step ithe mmaqEamet Ss tO réad a Sekaer 
input data and perform some preliminary calculations. 
This is accomplished by the main program, MAIN. 
ire col lowing sae ito eet menewLiput ; Ceteriinea 
from a Simple radial equilibrium design, which must be 
Submitted: 
(1) workang flumeeespeeitic gas constant sana 


specific heat ratio 
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MAIN: Overall 


no more 


read in data set - 
aata 





process some input data 
and print out selected values 





Li=1 


L2=1 
L3=0 


iterate to find a set of 
streamline positions and V; 
for plane I which satisfies 
continuity 





have 
~ yes treamlines ch ae 
L1=Li+1 streamlines changed 
much? 
no 









calculate and vrint 
velocity triangle data 





Figure 3-2 
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(2) turbine mass flow rate and rpm 
(3) total temperatures and pressures at all 
three planes 
(4) Shab and Meisecradamrrat Ole nae one and three 
(5) stator and rotor widths 
(6) coefficients of terms in the simple radial 
equilibrium relation for tangential velocity 
at each plane 
and (7) a reference velocity and radius used to 
fine the simple radial equilibrium axial 
velocity at each plane. 
In addition, the number of streamlines desired and 
certain test values are read in. 
mie foldewingeecaveulations ef a preliminary nature 
are performed: 
(1) test to ensure the overall flare is not 
zero and calculate flare angles 
(2) calculate the basic parameters used to 
apply the flare correction to the axial 
Velocity oi paracrapn 5.5.9) 
(3) calculate dimensional test values from 
Towle. (Pal agiscindg Oo... .4 3) 
Da) wealicitate blade terghies 
(5S) calculate certain other often used combina- 
tions of input which remain constant for 


the data Set under consideration 


aa 





selected parameters are printed out prior to further 
execution of the program. 
Bee.) Calculating a Streamline Location 

mie Neart OF the whole program is the calculation of 
streamlines; the most basic step in the streamline 
generation is the calculation of the "next" streamline 
location, RNEXT, from the position of an adjacent stream-~ 
line, R(I,J), and values VXZ and VTZ, representing average 
values of the axial and tangential velocities in the 
channel between RNEXT and R(I,J). This calculation is 
performed by subprogram RNEXT(I), whose flow chart appears 
ter agure. 3-3. 

The total velocity, V, is determined from the axial 
and tangential components, and with equation (2.13) 
determines the static to total temperature ratio. The 
temperature ratio and the equation of state for a perfect 
gas, equation (2.10), are used in another compressible 
BMow relation, equation (2.17), to find a density RHO. 
Then the continuity equation (2.35) can be used to find 
a flow area, AREA. From 


RNEXT = / AREA ROE, )) @ (3.1) 


the location of the next radially outward streamline can 


be found. 


oy 








Subprograin RNEXT 


_ DM 
~ RHO-VXZ 





Figure 3-3 
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Since RNEXT depends upon average values of axial 
mocerangential velocities, which are in turn functions 
amatus, 1teration 1S requmred. This iteration is 
accomplished by the subroutine STREAM shown in Figure 3-4, 

STREAM utilizes RNEXT as a guess for the next 
streamline and then calculates the average velocity 
values. The average values VXZ and VTZ are then used 
to calculate a new RNEXT andes on. When finally the 
difference between any RNEXT and the previous RNEXT is 
small, STREAM returns the values R(I,J+l) = RNEXT, 
eee stl ieand Vitt,J+tl) to the main program. The index 
L3 counts the number of times RNEXT is calculated during 
the calculation of a set of JJ streamlines at any plane. 

Pieadagitrton to famcdang an R(1,J+l1), STREAM alse 
applies the flare and actuator disk corrections to the 
unmodified (Simple radial equilibrium) axial velocity. 
The corrections are accomplished by multiplying the 
unmodified axial velocity, calculated in subprogram FVX(TI) 
iereacranph 3.3.6), by the corrections for flare and 
actuator disk effects which are calculated by subprograms 
Pigerig(t) and DISK(I), respectively. (See paragraphs 
Beoeo and 3.3.10) Since the very first time a set of 
streamlines is Calculated at A plane, DISK is undefined, 
only the flare correction can be utilized in the first 


Gee) 1OOp. 
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N3=0 
RR = best guess of R(I,J+1) 
RGUESS=RR 
VX(I,J+1) = FVX(RR) 
VI(I,J+1) = FVT(RR) 





Correct VX(I,J+1) for flare and disk 


VXZ = 4(VX(I,J) + VX(I,J+1)) 
VIZ = $(VI(I,d) + VI(I,J+1)) 


Bie dca) SUMS On MTL Ara OARISHGIINA 
RR = R(I,J+1) 

VX(1I,J+1) = FVX(RR) 

VT(I,J+1) = FVT(RR) 






Hi 






= VX(I,d+1 )“FLARE(I) 


VX(I,J+1) = VX(I,J+1 )*DISK(TI) 


Figure 3-h 








RGUZSS=2(I,J+1) 
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foro «6 «Calculating ee and Vo (Simple Radial Equilibrium) 
In a simple radial equilibrium preliminary design, 
ims standard procedure to pick a functional form for 
the radial variation of tangetial velocity. Then, the 
mre Variation of the axial velocity is determined by 
the simple radial equilibrium equation (2.23). In the 
present case, the tangential velocity can have the form 
of equation (2.27) where the coefficients aot Ayr ave and 
b are chosen by the designer. The simple radial equili- 
brium equation has been integrated for this form in 
Pepenaix Bb, resulting in equation (2.30). The function 
Submsegrams FVX(I) and FVT(I) calculate the simple radial 
equilibrium values of axial velocity and tangential 
velocity, respectively, given an argument RR and input 
parameters A2(I), Al (I), AO(L); B(L), VAO( 1), ance (2) 
3.3.7 Generating a Feasible Streamline Set at a Plane 
The main program calculates a set of streamline 
Wecations at a plane starting with the hub and using 
STREAM to step outward to the tip. Because the axial 
velocity has been modified, however, the average axial 
velocity is different from the simple radial equilibrium 
design value, and the last streamline does not fall on 
miemetxeaq tip radius. To ,correct, this, a relerence 
Peal velocity, VXR(I), used by FVX(2), is adjusted and 


the streamline set is recalculated. 
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Mice Teeration is Continued until the last streamline 
Mecatt@n 1S close enough to the desired tip radius. A 
EelowecNnart of the section Of "MAIN that generates a set of 
streamlines that fit the geometry is shown IME Lounes Sor 
mie Variables S and Sl do not play a direct role in the 
streamline generation and will be explained in paragraph 
(3.3.11). The adjustment of the reference axial velocity 
is explained in the next paragraph. The index L2 counts 
the number of iterations made at a plane to converge to 
jae proper R({i,JJ). At the users option (by specifying 
KK=1 in the input data) L2 and other parameters may be 
printed out after each trial set of streamlines. 

Seeseo, Ad }uUSting VXR{1) 

Initially, VXR(I) is set to VxX@(I), the value 
determined from the simple radial equilibrium cesign. 
Whenever the location of the last streamline, R(I,JJ), 
messes the tip radius, RT(I), VXR(I) is adjusted to 
increase or decrease the average axial velocity before 
recalculating a new set of streamline locations, in order 
tommare R(Il,JJ) converge to RY{1). 

To find an expression for a correction to the 
reference axial velocity in terms of the error, 


DR = R(I,JJ) - RT(I), an expression for the flow 


= V_°A (3752) 


geo ee 
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MAIN: Iteration for Streamlines at a Plane I 


Calculate VX(I,1) 
and VI(I,1) 


L2=L2+1 S=0 
Modify VXR(I) DO 15 J=1,dd = 


PaLrs 11, 12,13,1 





| 
| 
CALL STREAY | 
(calculate new R(I,J+1) | 
| 
| 
| 








Is 
no RT(I)-R(I,JJ)| 
less then 


Peotze 






yes 


Figure 3-5 


She 








can be differentiated, giving 


mM — 
d (7) = Adv, + Vda. (3.3) 
But we desire ad (m/p) = (0) euigtol 7\ Foe ip USUI ces Pear aac 
whence 


2V RT (I) *DR 
dv, = aC (3.4) 
Puu(T) “eR (1, 1)°° 
From experience, to make the convergence more rapid, the 
following term was added to the VXR(I) by MAIN: 
= 2 VX(I,JJ)_DR 
avy ~ RP(I) - R(I,1) (3.5) 
feos COrrecting Axial@Velocity for Flare 
The simple radial equilibrium axial velocity 
calculated by subprogram FVX(I) is corrected by multiplying 
it by a non-dimensional factor computed by subprogram 
FLARE(I). The basic relation involved is equation (2.33), 


rewritten as 


FLARE = cos*¢ = —————— , (3.6) 


where 


RR_7_REP egh, 


I 


tan $ 
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beets the radial location of the point of interest, D is 
a parameter depending upon the plane under consideration, 
and REF is a geometric parameter. For additional details 
see Appendix F. FLARE is always positive and less than 
or equal to one. 
ae LY )6Correcting Axial Velocity for Actuator Disk 
Effects 

Miesaceuatou, disk correction 1S applied after the 
HectboeCOrrection and 1S also a multiplicative non-dimen-— 
Sional factor, calculated by subprogram DISK(I). The 
Domwetons used in DISK{I) are forms of equations (B.1), 
Oe veend (h.2) for planes one, two, and three, 
respectively. DISK is always positive and may be 
greater than one. 
eeo.ll Testing for Overall Streamline Convergence 

The corrected axial velocity at a streamline depends 
upon its location and the value of axial velocity on the 
same streamline at planes up and downstream. Therefore, 
in general it is necessary to recalculate the streamlines 
at each plane several times, the streamlines moving less 
with each iteration. The streamline poSitions are 
calculated in order at planes one, two, and three at 
Meast twice, Each time the streamlines for plane Chree 
are recalculated, the following sum is accumulated: 

JJ 


S = BNE rap 
J=1 


Rts 7.) 


this iteration _ last terati2en 


(330)) 
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Berer €O Figures 3-2 and 3-5. <f S is small enough, the 
streamline positions have not changed much from the 
previous iteration. This means that the set of corrected 
femal velocities, VX(I,J), and streamline positions, 
R(I,J), have converged to mutually consistant values, 

and iteration can be terminated. 

Peele «=Velocity Triangle Calculations 

Once the streamlines have been located at each 
principal plane, a velocity triangle can be constructed 
pemeacd pointa. Pl, P2, and P3 in Figure 3-1 are the 
three points for one particular streamline. JJ triangles 
are constructed at each plane. 

The fundamental input variables to the velocity 
metaniglie calculations are the streamline location R(2l,J), 
cies ana VX(I,Jd), which were calculated concurrenely 
in the iteration process. The following derived 


quantities were calculated by MAIN at each plane: 


U = 27 * RPM * R(I,J)/60 (3.9) 

v= YWT(I,5)2 + VX(I,d)2 oo) 
w= YVX(1,J)2 +[U - VT(I,J) ]? (3a) 
a = arctan (Soret) (Bel ) 
g = arctan (So) 2) ( Sraalee)) 


sal 





M = from equation (2.16) 


= 

I 

ce 
<I 


creuk (32.04) 


In addition, the following quantities were calculated for 


each streamline: 


V1 _ V(plane 2) 


V2 ~ V(plane 1) (3.15) 


W2 ._ W(plane 3) 


Wl  W(plane 2) (3.16) 


stator turning angle = a(plane 2) - a(plane 1) (3.17) 


rotor turning angle = 8(plane 3) - B(plane 2) (3.18) 


To compare the above results to those attainable by 
Simple radial equilibrium, all of the above velocity 
triangle calculations are repeated at the same locations, 
but with the unmodified simple radial equilibrium axial 
Vemrocity. 

Figure 3-6 is a flow chart describing the calculation 
and printing of the above results. 

PeeelS Iteration Test Values 

For flexibility, the user supplies the test values 

used to terminate the three iteration loops. The three 


loops are (1) overall streamline convergence, 
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MAIN: Calculate and Print Velocity Triangles 


bo LO -K=7,2 


WRITE heading 


\ 


DO. 35) 1=1,3 = ee 


WRITE he vo mmetes | 


TN 
.f 
\ 
= 


DOD 35 J=1,dd SSS 
Gao) = Fvx(1) <> 


, 
al 
a 

vr 
ot 







“i? 


WRITE data row J for plane I 





j = 
3 
es 


Calculate velocity triangle data for J 






WRITE heading 


DO Ova d ee ee, eee 


WRITS row J of interplane data 


ho 


Figure 3-6 
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me convergence of the JJ'th streamline at a plane to 
RT(I), and (3) the location of the next streamline 

R(I,Jj + 1). These three loops are counted and labeled 
Peete Variables 11, L2, andes, respectively. Test 
values TEST1, TEST2, and TEST3 are calculated from input 


parameters Tl, T2, and T3 and the following relations: 


teoti = Tl (RT Oy — R(2,1)) (3.19) 
esi 2 ee ee RT (2) ce RX 2? 1) ) (3220) 
test s= 9s Re (2) = R(2,1))/IM1 (3.2200) 


An iteration loop is terminated normally when the error 


term becomes less than the appropriate test value. 


3.4 Convergence Checks 

If for some reason the program ever begins to 
calculate diverging values of the error term in any of 
the three iteration loops, the program will automatically 
terminate the calculations on the current data set. This 
is accomplished, as depicted in Figure 3-7, by exiting 1£ 


the error term has increased at least twice in that loop. 


Ceemer>rogram Verification 
In this thesis a model to improve Simple radial 
equilibrium design by applying flare and actuator disk 


corrections to axial velocity is described and a 
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Test for loon Divergence 








Enter 
Mteravion 
Loop 







N=0 






Sone large value 


Regular loop 
- 


miscellaneous calculations; 
Lind new ERR 


| 
| 
| 
| 
| 
| 
| 
| 
| 





EBRROLD=ERR Fé 


no 
Wale. error 


message 


Figure 3-7 


6 








computer program to implement the model is presented. 

The computer program has been verified to the extent that 
there is confidence that it accurately reflects the model. 
Full evaluation of the model, including comparison to 
more complicated, exact design methods, must be left for 
future efforts. 

Although the iterative nature of the streamline 
calculations makes hand calculation of the whole procedure 
very burdensome, it is easy to verify a solution by hand. 
The outputs of subprograms FLARE(I), DISK(I), FVX(I), 
and FVT(I) were all verified by hand, given a computer 
calculated streamline position R(I,J + 1). With these 
values and R(I,J), RNEXT(I) and STREAM were verified by 
Nen@ecalculating R(I,J + 1). The computer calculations 
of flare geometry, velocity triangle related output, ana 
other algebraic relationships, all agreed with hand 
ealculation. 

Since each input data set must be a complete simple 
radial equilibrium preliminary design, the number of 
independent geometries, working fluids, tangential 
wellocity distributions, etc., tested was not very large: 
Mevertheless, it is felt the verification of the program 
was sufficient to ensure accuracy for feasible designs. 


Results are shown for three sample cases in Chapter 4. 
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A aaRESULTS 


4.1 General 

This chapter presents the results of applying the 
computer program to three sample simple radial equilibrium 
designs. The computer output listing is included and 
plots of representative variables are provided for easy 
comparison between the simple radial equilibrium input 


design and the modified design. 


4.2 Free Vortex Design 

Table 4-1 shows the characteristics of an axial 
inflow and outflow free vortex stage designed as the last 
stage of a double flow turbine in a geothermal power 
peiefes tion.” A sketch of the geometry appears in 
Figure 4-2. Notice the design has a constant hub radius, 
but has large flare in the shroud. The working fluid is 
Freon 21. Eleven streamlines are called for. 

The computer printed results are listed in Tables 
4-3, 4-4, and 4-5. It 1S seen that the large changes in 
axial velocity (about 25% at the hub) can result in a Six 
degree change in the relative inflow angle to the rotor. 
Less pronounced, but readily noticed changes are shown 


in other flow angles, velocity ratios, and Mach numbers. 
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ite igure 4-6, the axial velocity, stator outlet 
aagle, ©Ootor fluid turning angle, and relative velocity 
ratio are plotted for both modified and original simple 
radial equilibrium designs. It may be Significant se 
note that the changes at the hub (where ¢ = 0°) may be 
considered to be "indirect" in that they result from the 
increased average axial velocity required to obtain the 


same mass flow. 


4.3 Constant Reaction Design 

The input parameters for another example having the 
Same annular geometry as above but with a different 
assumed tangential velocity distribution are shown in 
Table 4-7. Axial inflow to the stage is maintained, 
but the constant reaction type simple equilibrium 
relations (see Appendix D) are used at the rotor leading 
and trailing edges. 

The computer printed results are shown in Tables 
fageeeo-9, and 4-10, and plots of “axial velocity and rotor 
relative flow angles appear in Figure 4-1l. From Figure 
tai) the Penn eation to axial velocity doesn't seem as 
significant as in the first example; the influence on the 
rotor flow angles is still significant, however 


(AB = 6° at the hub). 
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4.4 Free Vortex, High Mach Number Design 

To demonstrate the adequacy of the computer program 
to calculate in regimes where M 2 1.0, the preliminary 
design summarized by Table 4-12 was run. In this case, 
the flare 1S symmetrically distributed between hub and 
shroud, as 1S shown in Figure 4-13. | 

The detailed results are listed in Tables 4-14, 
a>, and 4-16. The axial velocity and Mach number at 
the rotor inlet (plane 2) are plotted in Figure 4-17 


for comparison with the simple radial equilibrium values. 
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De eONGhl SONS 


5.1 General 

A model for improving a simple radial equilibrium 
design by modifying the axial velOoey has been ciscussed. 
The computer program presented here provides a means to 
quickly and inexpensively incorporate the model into 
preliminary design procedure. 

The improved preliminary design model is basically 
a perturbation of the simple radial equilibrium design 
procedure, and comparisons of the two procedures can be 
made easily from the computer program output. It would 
be useful to compare the improved model results to designs 
produced by more rigorous and complicated analytical 
procedures as well, to see where the balance of complexity 
versus accuracy lies. Since the program described here 
is only a "design tool" and was not derived analytically, 
comparisons to more exact procedures are really needed to 
eeeve full respectability to the results. 

The changes made as a result of applying the model 
tO a Simple radial equilibrium preliminary design are of 
a magnitude which suggests that they may have measurable 
effects upon turbine performance. Further study is 


needed to quantify these effects in terms of efficiencies. 
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5.2 Improving the Model 

Aside from possible improvements in the generation 
Seeene Correction to the simple radial equilibrium axial 
velocity, other changes in the program to extend its 
usefulness may be desirable. An example is extending the 
range of the program to include more than one stage, 


with two conical sections per stage. 


wo. Applications 

The most direct application for the program is to 
refine the blade shapes of an axial turbine in preliminary 
design. However, less specific applications to research 
type problems may be fruitful. Used in conjunction with 
a sufficiently general simple radial equilibrium 
preliminary axial turbine design program, the program 
presented in this thesis could be used to generate a 
series of modified designs to sie the ramifications 
of flare and actuator disk effects. Perhaps it is 
possible to find a simple radial equilibrium tangential 
velocity distribution which will be less effected by the 


Meer tication, Or will “compensate” tor the etfects. 
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APPENDIX A 


Variable List 

The following is an alphabetical list of variables 
which appear in the computer main program and subprograms. 

The units are the same as those listed in the 
Notation section. The same variable names are used to 
represent identical quantities in both the main and sub- 
programs. In addition, the subscripts are used in 
consistant manner throughout. An intermediate variable 
mmone used locally in the program to simplify arithmetie 
expressions. Zero is written @. Unless otherwise 
indicated, variable naming follows the implicit REAL/ 


MieGeR typing convention. 


AREA annulus area between two adjacent 
streamlines 
Ag (1) constant term in the expression for 


tangential velocity 


ja Gal) ; coefficient of r in the expression 
for tangential velocity 


A2 (TI) coefficient of r* in the expression 


for tangential velocity 


B(I) coefficient of 1/r in the expression 
for tangential velocity 


ik INtermediate Wace PuOpeteyr colotant 
used by RNEAT 


CZ Intermeatate variable 


OS intermediate Variable 


20 








D 
DALPHA 
DBETA 
Ph{t,J) 


Beok (1) 


DM 


DR 


DROLD 


DRZ 


DRZOLD 
ee, 1) 
Fac yc ) 
Hea, 3 ) 
F(I,4) 
FF (I,J) 


PLARE (1) 


FLOW 


Ey? (1) 


BVE (LT) 


GAMMA 


intermediate variable 

Stator tuning angle 

rotor turning angle 

actuator disk correction factor 
function subprogram which calculates 
Che snen=cd MemsaOiN a eeaetel ah erence clc 


correction factor 


mass flow per channel between 
streamlines 


difference between the calculated 
and actual tip radius 


value of DR in last iteration 


difference between R(I,J+1) and RGUESS 
Galeutated by SPREAD 


value of DRZ in last iteration 

total absolute gas velocity 

total relative gas velocity 

alpha; absolute flow angle 

beta; relative flow angle 

flare correction. tacrer 

function subprogram which calculates 
the non-dimensional flare correction 
Tacetou 

total turbine mass flow 

function subpreqreamevimel calculates 
tangential velocity as a function of 
radius 

function subprogram which calculates 
axial velocity as a simple radial 
equrlibrium 2unetlOnmcremadUS 


g, = 32.1740 £t-lbm/1bf/sec* 


y = ratio of specific heats; input 
constant 
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HR 


HS 


16)8) 


IREAD 


IWRITE 


J 


JJ 


JM1L 


KK 


i 


L2 


1 ES: 


MV 


mean rotor height 
mean stator height 


subscript which identifies a transverse 
Station along the turbine axis 


1 stator leading edge 
2 rotor leading edge 


I 
if 
me 3 rotor trailing edge 


No fl 


user determined input data identification 
number 


input device channel number 
Output device channel number 


streamline index; J = 1,2,3, ..., Jd, 
starting at the hub 


quantity of streamlines considered 


quantity of flow channels considered 


key variable; K 1 mocified simple 
racial equi 


brium data 
Ko = ee okt, data 


key variable; KK # 1 program diagnostics 
not printed 


KK = 1 diagnostics 
printed 


loop count; number of times all the 
streamlines for the turbine are 
calculated 


loop count; number of times a set of 
streamlines must be calculated at a 
station to satisfy the geometry 

loop count; the total number of times, 
Quming the calcilatlonpOt wd Set Ou 
streamlines at a plane, that the 
subprogram RNEXT is called 

key used in labeling output data 


absolute mach number (REAL variable) 


a2 





MW 
PUB 
PML rEP 


Nl 


N2 


N3 


el 
a2) 
Ear) 
REG 
RGAS 
RGUESS 
RHO 


RNEXT (I) 


RPM 


RR 


Rt) 


age) 


SOLD 
Sl 


A Uierpevg bab 


relative mach number (REAL variakle) 
hub flare angle 
tip flare angle 


test variable for Li loop divergence 
check 


test variable for L2 loop divergence 
check 


test variable for L3 loop divergence 
check 


Sails so) 

stagnation pressure 

radial location of a streamline 
reference radius for flare correction 
specific gas constant; input data 
temporary storage for a trial R(I,J+1) 
static density 


function subprogram which calculates 
a tCrialavValMeror hn, wrt) 


turbine rpm; input data 


intermediate variable used as the 
argument fOr FUARE GE. Vet), EVE as 


np bacwus 


input reference radius at which VX(I,J)= 
VXp (I) 


Sum When accumilatesetme Gailrerence 
between streamline positions on the 
current iteration (Ll) and the 
previous (L1-1) 

value of S from previous iteration 


teria meee 


test value to terminate loop Ll 


els 





wei 2 
eo! 3 


TEST4 


TGNT 
TR 
we 1.) 
De 
T2 


In 


VRATIO 
Wee Id) 


ee 


wet, J) 


VXR(I) 


Meacnr (1 J) 


VXZ 


VXP (I) 


WE 


WR 
WRATIO 
WS 

Z 


7 Oe ary 6 


test value to terminate loop L2 
test value to terminate loop L3 


CeSte- Valuecieo ensure thera). 
velocity is positive 


tangent of the flare angle §@ 

GatlOwon Statice gke COLaImeECmMbperarurc 
stagnation temperature 

input percent used to calculate TESTI 
input percent used to calculate TEST2 
input percent used to calculate TEST3 
blade speed 

ae: 

tangential velocity 


average value of ‘tangential velocity 
between two streamlines 


axialtveloci ey 


reference value of axial velocity; 
input data 


simple radial equilibrium axial velocity 


average value of axial velocity between 
two streamlines 


input reference value of axial velocity 


entrance wicath dimercrven sea to Line 
theswtlare correction teactor 


YOtOr wile cm 

VALET 

Stator wlan 
intermediate variable 


intermediate variables 
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APPRENDI. 


Program Listing 

This appendix contains the computer program listing 
of the MAIN program, subroutine STREAM, and function 
Subprograms RNEXT, FLARE, DISK, EVX;, anG@ FV Eee eter Ito 
the table below. The purpose of each program segment 
is indicated; flow charts for MAIN, STREAM, and RNEXT 


appear in Chapter 3. 


paragraph seoment page 
Bet MAIN 35 
Bie STREAM 10K 
B38 functions 106 


B.1 MAIN 

The MAIN program reads the input data, governs the 
iteration loops to generate the streamlines, calculates 
mae Velocity triangles, and prints ae the results. The 


listing begins on the next page. 
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BeZ2 Subroutine STREAM 


Sinem calculates an ‘x. iterating on assumed 


diy tl 
memenes for V. and V, until xr. . and the average 
x S) dee ato 
properties based on r. . andr. . are consistent. 
Ly 3 i ee 
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Pee Function Subprograms 
The listing for each of the five function subpro- 


gramS appears on one of the following pages. 


Bunce Lon Page Description 
RNEXT 107 given a value of rs : and the 
f 


properties ex Vor RNEXT 
Caleculaves san tlic lee 
veyed fart 


FLARE 108 given r. FLARE calculates a 


TF gill’ 
non-dimensional factor which, 
multiplied by Ves in STREAM, 
gives a Ve modified for FLARE 

DISK 109 given appropriate values of 
en! DISK -GalecuLtates a None 
dimensional factor used to 
modify axial velocity for 
actuator disk effects 

aval 110 given ee ae PVT ee abewilates 2 Vo 
by a proposed simple quadratic 
relation | 

FVX Eas given raat FVX calculates a V, 


using a Simple radial equili- 


brium relation 
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AE EN Dix ae 


Beecitic Instructions For Use of the Program 
C.l1 General 


This appendix is intended to provide detailed 
instructions for the use of the computer program 
developed in this thesis. Preparation of the input 
deck and interpretation of the output will be covered. 

The program was run on the Interdata 70 computer 
at the Joint Mechanical and Civil Engineering Computer 
Facility at MIT, but should be easily adaptable to any 
BoRURAN IV capable system. Keterences 7 "and 67 served 
as the FORTRAN IV language guide. All READ and WRITE 
instructions are written with dummy variables IREAD and 
IWRITE, respectively. Therefore, the only two FORTRAN 
statements which need to be changed in order to run on 
another system using different input/output’ channel 
numbers are the two assigning integer ates to IREAD 
and IWRITE. The control cards, listed below, which may 


be system dependent, are marked by an asterisk. 


C22 Loading Sequence 
Listed below, in the sequence that they should be 
loaded, are the contituents of the program: 
mL) BOD ease 


*(2) Language and other compilation control 
cards 


bby 





(3) MAIN program 


(4) Subprograms 


a) 
b) 
Cc) 
d) 
e) 


t) 


m5) =Chxecute card 


(in any order) 


subroutine STREAM 


FUumGtlen 
hunict en 
LUnee LON 
£Unieta on 


1B 0a Yoxenitay el 


subprogram 
Subprogram 
subprogram 
subprogram 


Subprogram 


RNEAT 


FLARE 


DISK 


Ey. 


- 


tvs 


*(6) Common length designation and the 


execution control cards 


(7) Data cards 


*(3)° JOb terminations cane 


Gms cob Control Cards* 


The program segments marked by an asterisk may 


vary in content and format from system to system or with 


time, and must be determined from an appropriate user's 


manual. it should be nored that tie amount oo. bilan 


common utilized by the program depends upon the dimension 


mercome of the variables, which relates to. ene numbersor 


SEucamlines Conslaerec. 


C.4 Data Cards 


The data deck is composed of any number of separate 


data sets stacked directly together. 


Each set is 


independent of every other set and represents a separate 


simple radial equilibrium turbine design to be processed 


by Ene program. 


Nes: 





Fach data set consists of ten Separate cards which 


must be present and in the proper sequence. A specific 


femmat 1S given for 
format (data fields 
the computer system 

The context of 
meacagraph C.4.1. A 
note in paragraph C 


Gea.t Format Table 


Columns Format 
CARD ONE 

4 r4 

She 14 
9-12 I4 
CARD TWO 

lO F10.4 
ie 2.0) PLO. 4 


eaRD THREE 


7-9 |S 
io-18 Hoes 
io] Be S s: 
Bo 36 SAS) 
B7=45 | io ee 
46-54 Bo. 3 
52-63 Bos 
p47 2 ou 


the data on each card, although free 
separated by commas) may be used if 
has that feature. 

each card iS given below in 


Superscript refers to an explanatory 


ee. 

fneornat Lon Dine s 
number of Seman nos- none 
diagnostics print key~ none 
user data set ame none 
gamma (ratio of specific none 

heats) 

specific gas constant ft Ibf/ lbn7ee 
turbine rpm rpm 
total mass flow rate lbm/sec 
total temperature, plane 1 oR 
total pressure, plane l psi 
total temperature, plane 2 OAR 
total pressure, plane 2 psi 
total temperature, plane 3 YR 
total pressure, plane 3 psi 


114 





CARD FOUR 
= 10 
LiSZ0 

Za 3.0 
B40 
e150 
Sao 0 
CERD FIVE 
TO 

= 20 

Zale 3.0 
31-40 
GARD SIX 


J=20 


iM 20 


Zi 30 


31-40 


ee 20 


Dee 6 0 


CARD SEVEN 


0 


tek 20 


Zales 


31-40 


Bao. 
F102 
| BO 
Faron 
ie OP 


EEO s 


gd JO 


EAEOR 


Ear, 


FLO: 


ie IOs 


BO 


BO 


lade) 


F10. 


EO 


jel 8 
EEO. 
Per 


BO 


hub radius, plane 1 
tip radius, plane 1 
hub radius, plane 3 
tip radius, plane 3 
stator width 


FOOL WiloEn 


test percent i 
test percent 2 
ieg(eie (Ceneciennie 3 


test value 4 


SRE reference radius, 
plane l 


SRE reference velocity, 
plane l 


SRE reference radius, 
plane 2 


SRE reference velocity, 
plane 2 


SRE reference radius, 
plane 3 


SRE reference velocity, 
plane 3 


8 bet jelsuye ee 
AG FOG wlancw 


A, for plane 1 


A. for plane l 


AS 


Be 
ine 
Zt 
ft 
ie 


IIe 


GGG Ie 


none 


none 


ce 


Ft 


ft/sec 


Ete 


tt) see 


iets 


ft/sec 





eaAnD EIGHT 
fsame aS card seven, but for plane 2) 
CARD NINE 
(Same as card seven, but for plane 3) 
Carn) TEN 
(The tenth card is blank) 
C.4.2 Notes 
1. The limit on the number of streamlines is only 
a practical one, related to machine computing speed and 
storage, and can be any reasonable number. The program 
was Originally run with eleven; any greater number will 
require a new DIMENSION statement for the main and all 
subprograms. Additional blank common Space may also 
have to be allocated. 
2. Tf this data bit=1, allvavamleblesdiaqnesceres 
will be printed out. See paragraph C.5.3, below. If not 
equal to unity, no diagnostics except error messages will 


be printed out. 


3. The user may use any three digit number to 
identify particular data Sets. “Some mumber must be 
Supplied. 


Wee ReLGT EO paragrepl ast o- 
Se. 5, ae Ay: Ayr are the coefficients of powers 
©f © in the tangential velocity relation. See 


bemaqraph 363.6. 
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Gee Output DUSCription 
C.5.1 General 


The output for each data set may be considered in 
three parts: 
Cem int rodute cony 
(2) diagnostic (optional) 
(3) velocity triangle data 


(1) modified simple radial equilibrium 
values 


(TAI wGCRE Values. 

See Introductory Output 

Input data and other selected calculated values, 
which do not change during the execution of a data set, 
gieeprtnted Out in the introductory block. The Symbols 
appearing have the same meaning as used in the program 
el@egiven an Appendix A. If the datawis anout lalicc, 
the Output units are the same. TESTI, TEST2, TESTS, 
REF, and WE have units of feet, and PHITIP and PHIHUB 
are in degrees. 
G5.-5 Diagnostic Output 

If input variable KK=1l1, certain lines of information 
are printed out during execution of the program iterations. 
Sample output of this type is shown at the end of this 
appendix. Each line contains the current values of 
certain important variables of the program at a certain 
step in the calculations. This information was used in 


the development of the program and should not be 


belly 





required routinely. However, it can increase confidence 
in the results, and may be useful if unexpected problems 
@eeur during operation in the future. @Additionalecactails 
on the meaning of the loop counters can be learned from 
Eme ttow charts presented in Chapter 3. 
Gao-4¢ Velocity Triangle Output 

A group of results is printed for each of two cases: 

(1) The streamline numbers and locations, and 
velocity results are printed for a simple radial equili- 
brium preliminary design which has been modified by the 
flare and actuator disk corrections determined by the 
program. Results pertaining to each of the planes is 
printed, followed by a block of results involving values 
from more than one plane. 

(2) For comparison purposes, values calculated 
at the same streamlines as above, but with no modifica- 
mon for flare and actuator disk ame printed aeareie 
same format. 

The meanings of the variables are the same as 
festeadsin Appendix A, Or are self-explanatory. iste: J 
and 2 appearing in some column labels werfer to rocor 
inlet and outlet triangle, respectively, following the 
usual convention. Distances are in feet, velocities in 
ft/sec, and angles in degrees. The value VXR(I) printed 
for each plane is the final value for the constant term 
in the FVX relation for axial velocity. It was changed 


from VX@(I) to satisfy continuity. 
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C.5.5 Error Messages 

The input data to this program is assumed to be 
from a feasible preliminary design. Therefore, many 
mypes Of execution errors; stich as fromeene appearance 
of negative radii, should not occur. If some inconsis- 
tant or infeasible data should be used, however, measures 
have been taken to minimize the chances of hanging the 
Beegram up in an iterationM@ioop. The ceetails of the 
testing procedure to accomplish this are found in 
paragraph 3.4. Each iterative loop is provided with a 
means of exit if the error term increases for at least 
two successive iterations. Another check is made to 
determine if the axial velocity is less than 1.0 ft/sec 
(an arbitrary, small positive value, to avoid the 
possibility of Sse sone: bY -ZEre)- 

If any of the conditions listed in the next para- 
graph are discovered, an error message will be printed 
out regardless of whether or not the user has specified 
=e printing of the diagnostic dara. 

The following is a list Of error conditions tested 
for, and their corresponding output messages: 

(1) Condition: the axialmeelocity computed way 

the function subprogram FVX is less than 
1,0 Et/sec 
Error message: Ll = xxx L2 = xxx L3 = Xxx 
i = XXX d= ee”. VA (lad ese 


THAN eC Ew SEG 


Ae 





(2) 


=) 


ce, 


Action: execution continues 

Condition: the change, S, in the I=3 plane 
streamline positions is nor 
decreasing on Sieeeseiue iterations 

Error message: ITERATION TERMINATED DUE TO 

INCREASING VALUES OF § 
Action: execution of current data set is 
terminated 

Condition: the annulus geometry error, DR, 
is not decreasing on successive 
iterations 

PieGor eines Sage: TTERATION LeEkRMINATED DUE “TO 


INCREASING Ahr S TOF DR 


TL =" Se 2 ee se a= 
Ae Torn: execution of Current cata Set as 
terminated 


Condition: the -change,, a7. in the Calculate won 
Of R(l,Jt)) by sueueutine STREAM 
is not decreasing on successive 
iterations 
Error message: |) ITERATION TIE RMINATED DUE 10 
DIVERGING POSITION OF STREAM- 
BLING See I = xxx J = XXX 
LL = xxx G26 = xx 3S) = Sexe 
Action: execution of current data set 1s 


terminated 


dh 9) 





C.6 Additional Options 

In some applications it may be desired to find the 
Simple radial equilibrium streamlines, or to consider 
the effects of the flare or actuator disk correction 
alone. To accomplish this, the executable statements 
in one or both of the subprograms FLARE(I) and DISK(TI) 
may be replaced by the statement FLARE = 1.0 or 


Drsk = 1.0, aS appropriate. 


lak 





weir res 


asd) iy Sag 


rh b “ey * 


Oe eden rae (a. 


tol rege 


& dee ee 
iis iat a Se! 
Gt 121° 
ct 2" 
Cauley cal 
as Fy? 
Vole? 


ban ss d*tel 
*stwie%tad 
"4YG4O3 ob 
ls€ ’e 3% ct bb 
1 z(°M el] de 
sunu“Le*ted 
l =a(TCeldy 
Z =(fC ella 
*Jo4Gs°% tok 
Di ae ae Ga Ee 
gyodlucoront 
7 2( PO ey de 
t i hee aa 0 ae 
L =F elle 
CVV e2a° wos 
if 2 (relia 
l e(TMetia 
b-E/es° fol 
ak z¢rUelly 
> 2zIiTlMel day 
Lb F(F Mel )g 


=C1)eaxa 
=(T)exKa 
=CLdI-AaIXA 
sO] aka 
bleivet = 
=H=C1)exKA 
bTvnes 
wbhLeWd* 
=H] )cxa 
ale 
=CL 4A 
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liv aAeye 
GEV Beg 
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bev jevyeo 
LLGAe¢, 
=(1l)>xa 
T*evianey 
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Lindet 


ivyj*G¢o2 


23144° 
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APPENDEX 5 


Integration of the Simple Radial Equilibrium Equation 
A convenient form of the simple radial equilibrium 


equation 1s 
V 
See Oe AG 


Given explicit functions, Vg (re), it can be easily 
integrated to find a relation for axial velocity as a 
Bumetion of radius. 

Suppose the expression co the tangential velocity, 
Vo (x), can be written in the following epee 


= 2 
Vg (x) = ajr° + ayr + a, + 


i (Gs 


(Dn.2)) 


The bottom sign applies after a stator row, and the top 
sign applies after a rotor row. 
Then, proceding to evaluate the right hand side of 


equation (D.1), 


a _ > + 

qr (Vg) = 3a5X + Za,r a 

eae ee (D.3) 
esac 8 2 1 r 


123 





miaally, multiplying (D.2) by (Dv3), (be becomes 





av 
i : Sie 2238 2 2 
ue ae 3a. r°> + 5a,aor + (4a a, + 2a, eg 
fe 3 2 i an 
( aoa 7 3a,b) + (a. 3 2a, b)= = 2 (D.4) 


Equation (D.4) may be integrated directly from a 
reference radius, rie and axial velocity, Vow to an 


arbitrary radius r, and corresponding Vee This gives 


} fa 2 Paes: 2 4 4 5 3 : 3 
—_ —_— - — — - — — A 
7(V V ) ga.“ (xr r.') + $a, a,(r or) 


+ (2a a, ta 
O 


2 2 2 - a 
9 1 (ie r. ) + 3 (ajay = ad) (x r.) 


Na an ae 
+ (a. a Ber aaa =: ab (> i 


Rearranging, 
ane Dene Opies 9 ee tee nga 3 3 
Vi, ae a, (x re ) “a, a, (x re ) 
ae 2 pee 2 nee = = 
(4a a, +2a, ) (x rs ) 6 (aia, m ab) (x r.) 
peo encase) 2 2A ote ee fish 
O af 7 or xr r. 


The general expression, (D.5), may be simplified 


for two special cases in common use: 
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(1) Let a5) =a ae ays b 7 02 Then). 5) becemes 


a 2 ey 2 
Ve Vw 2a, (ie rs Bs 4a,b In(r/r,). 
(D6) 


This represents the constant reaction type design. ° 


fr) Let a, = a, = 0; aye b #~ 0. Then (D.5) becomes 
2 aes eee _ tL. i 
i Ne 2a. in(r/r, ) 7 Zab (= ae 


(Died) 


This represents a design in which the stator inlet angle, 


Oo 1S approximately the same at all radii.° 
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APPENDIX E 


Actuator Disk Superposition 
Approximate solutions of the actuator disk theory, 
discussed in Chapter 2, are, in the case of the 


H=olated disk: 


Ove) = Ga) 


_ 400 00 Tr 
ve = Re + a i (CE) 
Mmeotream Of the @isk (x < 0); end 
Cnn = WU). ey 
CE 2 nearer A ee 


downstream of the disk (x > 0). 

These equations pertain tovemeracral  variaelenwes 
axial velocity along a streamline. 

Suppose it iS desired to consider the effects of 
these actuator disk solutions at the three principal 
planes of a Single stage axial turbine. Two actuator 
disks are involved, corresponding to the stator row and 
the rotor row. The three planes of interest are the 
stator leading edge (plane 1), the rotor leading edge 
(ellane 2), and the rotor) trailing vedgemmplane 3). 

For simplicity, assume that the stator leading 
edge is influenced only by the stator actuator disk 


(for equal blade widths, the rotor disk is three times 


AG 





the distance away), and that the rotor trailing edge is 
influenced only by the rotor actuator disk. The rotor 
leading edge feels the effects of both the stator and 
rotor disks. Refer to Figure E-l. 

The relations (E.1) and (E.2) apply to the stage 
Helet and outlet planes respectively. “The effect mat 
the rotor leading edge, point S, can be obtained by 
superimposing the velocity corrections due to the two 
disks. The effect of the stator upon S is given by 
(F.1) with x = -WS/2 and h = HS; the effect of the 
rotor is given by (E.2) with x = +WR/2 and h = HR. 

In addition, the velocities at infinity, appearing 
in the two equations must be interpreted. Let CV ere 
denote one axial velocity a long way downstream of the 
motor Gdisk? similarly, let OV aes denote the axial 
velocity a long way upstream from the stator, and so 


om. then (8.1) and (E.2) become 


(View) = (Vee) 
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Now assume 


(VO). = (VU)... = Vo. = the simple radial equili- 
ee ee = brium Vel@eecity at vene 
rotor leading edge 
(plane 2) corrected for 


flare; 
eS ee = Vow = the axial velocity at plane ike 
Dees = Ws = the axial velocity at plane 3. 


Then finally, the corrected axial velocity, V at 


sae 


point S (plane 2) can be written 


tT? WS 
Fe) BS, 





ee ie x+ Vem ae ass _ Vien ae (- 
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(Es) 
V9 1s the axial velocity correctedmiern both flare and 


actuator disk effects. 
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APP Ea 


Calculation of the Flare Geometry Factors 

The purpose of this appendix is to derive the two 
parameters REF and WE used by the computer program to 
eOmpute the Elare angle o>. The computation of @ at an 
arbitrary point, 0, within the annulus is also demon- 
Strated. The stage inlet and outlet hub and shroud 
radii are given information. 

Referring to Figure Pel, Gnwene next page, it 2s 
seen that the point P is the intersection of the extended 
line segments forming the hub and shroud walls of the 
turbine. The line x = 0 represents the inlet plane of 
the turbine, x = Xo represents the outlet plane, and 
y = 0 represents the axis of the machine. 


The equations for the lines representing the shroud 


and hub walls, respectively, are written as 
y = —————_ x + R (Fw) 
CE ee 


Solving (F.1) and (F.2)) simuttancous hy sie. x andy, the 
coordinates of the point P, namely (-WE,REF), are found 


to be 


IEENe 





(X3,R, 4) 
Shroud 3? t3 

Oe Ry 4 ) O(%55Vo) 

P(-W, REP) d 
ON i ae | ee ieee. | ee 
(0,R4) 
a ae (35,3) 
—— wo 
a x 
Figure F-1 
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(R => Roe x 
WE = qo 1 WE ? 0) Ci 3) 
Eo tl heen 
= WE z 
REF = Rol - x, (Ria Ria): (F.4) 


The angle ¢ for some arbitrary point, Q, within the 


annulus, can be found from the expression 


Y, ~ REF 


a hE ee eee 
O 


where the coordinates of the point 0 are (Kor ya 


Le 
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